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Abstract: We summarize recent work on the consistent calculation of bubble-nucleation rates. Our 
approach is based on the notion of a real coarse-grained potential. The bubble-nucleation rate is 
calculated through an expansion around the semiclassical saddle point associated with tunnelling. We 
resolve outstanding problems related to the convexity of the potential, the double-counting of the 
effect of fluctuations and the inherent ultraviolet divergences. We determine the region of validity of 
the expansion around the saddle point. We find that this expansion fails near the spinodal line, and 
for weak or radiatively induced first-order phase transitions. We apply our method to the bound on 
the Higgs-boson mass from vacuum metastability and the electroweak phase transition. 
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Developing a consistent theory of bubble 
nuclcation is an important problem for par- 
ticle and statistical physics. Possible applica- 
tions include the description of most cosmologi- 
cal phase transitions, the questions of tunnelling 
in quantum mechanics and field theory, the prob- 
lem of vacuum stability and its phenomenological 
consequences, the (currently unsatisfactory) ex- 
planation of the experimentally measured bubble- 
nucleation rates for statistical systems, etc. The 
estimates of bubble-nucleation rates for first-order 
phase transitions are usually carried out within 
Langer's theory of homogeneous nucleation [1], 
applied to relativistic field theory in refs. [2] . The 
nucleation rate is exponentially suppressed by 
the action (free energy rescaled by the tempera- 
ture) of the critical bubble, a saddle point of the 
free energy of the system. Significant contribu- 
tions to the nucleation rate may arise from higher 
orders in a systematic expansion around this sad- 
dle point. The first correction has the form of a 
pre-exponential factor that involves fluctuation 
determinants around the critical bubble and the 
false vacuum. The evaluation of this factor is a 
difficult problem at the conceptual and technical 
level, as crucial issues associated with the convex- 
ity of the potential, the divergences of the fluc- 
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tuation determinants and the double-counting of 
the effect of fluctuations must be resolved. 



1. The method 

In a series of recent works [3] -[7], following the 
proposal of refs. [8], we developed a consistent 
approach, based on the effective average action 
Tk [9] that can be interpreted as a coarse-grained 
free energy. Fluctuations with characteristic mo- 
menta larger than a coarse-graining scale (q 2 ^ 
k 2 ) are integrated out and their effect is incorpo- 
rated in Tk- In the limit k — » 0, Y% becomes the 
effective action. The k dependence of Tk is de- 
scribed by an exact flow equation [10]. This can 
be translated into evolution equations for func- 
tions appearing in a derivative expansion of the 
action [11]. An approximation that is sufficient in 
many cases takes into account the effective aver- 
age potential Uk and a standard kinetic term and 
neglects higher derivative terms in the action. 
The bare theory is defined at some high scale A 
that can be identified with the ultraviolet cutoff. 
It is, however, more convenient to choose a start- 
ing scale fco below the temperature T, where the 
effective average action of a (3 + l)-dimensional 
theory at non-zero temperature can be described 
in terms of an effective three-dimensional action 
at zero temperature [12]. 
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In ref. [3] we computed the form of U k at 
scales k < ko for a theory of one scalar field by 
integrating its evolution equation [10, 11] 

A W ) = -ilr^W). (1.1) 

We considered an initial potential Uk with two 
minima separated by a barrier. Uk is real and 
non-convex for non-zero k, and approaches con- 
vexity only in the limit k — > 0. The nucleation 
rate must be computed for k larger than the scale 
kf for which the negative curvature at the top 
of the barrier becomes approximately equal to 
—k 2 [13]. For k < kf the form of the potential is 
affected by field configurations that interpolate 
between the two minima (similar to critical bub- 
bles). For k <; kf the typical length scale of a 
thick- wall critical bubble is <; l/k. 

We performed the calculation of the nucle- 
ation rate for a range of scales above and near 
kf. The nucleation rate is given by 



I = A k exp(-S k ), 



where 



A k = 
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Here <j>b(r) is the profile of the spherically sym- 
metric saddle point, Sk its action determined by 
using the potential Uk(<fi), and <j> = corresponds 
to the false vacuum. The prime in the fluctu- 
ation determinant around the saddle point de- 
notes that the 3 zero eigenvalues of the operator 
—d 2 + U' k \(f)b{r)), corresponding to displacements 
of the critical bubble, have been removed. E 
is the square root of the absolute value of the 
unique negative eigenvalue. The above form of 
Ak guarantees that only modes with character- 
istic momenta q 2 ^ k 2 contribute to the nucle- 
ation rate. (The effect of modes with q 2 <; k 2 is 
already incorporated in the form of Uk{4>)-) We 
emphasize that the use of an ultraviolet cutoff in 
eq. (1.2) that matches the infrared cutoff proce- 
dure in the derivation of eq. (1.1) is crucial for 
the consistency of our approach. In both cases, 
mass-like cutoffs have been used. 



2. Typical calculations 

The details of the numerical evaluation of the 
complicated determinants in eq. (1.2) arc given 
in ref. [3]. Fig. 1 exhibits the results of the cal- 
culation for an initial potential 
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that has a form typical for first-order phase tran- 
sitions in four-dimensional field theories at high 
temperature. Through a shift <f> — > 4> + c the cu- 
bic term can be eliminated in favour of a term 
linear in <j>. Therefore, eq. (2.1) also describes 
statistical systems of the Ising universality class 
in the presence of an external magnetic field. 
The first row corresponds to parameters m| = 
-0.0433 kl 7fco = -0.0634 k 3 /2 , X ko = 0.1 ko- 
The first plot shows the evolution of the potential 
Uk{<f>) as the scale k is lowered. (We always shift 
the metastable vacuum to <f> = 0.) The solid 
line corresponds to k/ko = 0.513 while the line 
with longest dashes (that has the smallest bar- 
rier height) corresponds to kf/k — 0.223. At 
the scale kf the negative curvature at the top 
of the barrier is slightly larger than — k 2 f and we 
stop the evolution. The potential and the field 
have been normalized with respect to kf. As k 
is lowered from fco to kf, the absolute minimum 
of the potential settles at a non-zero value of <f>, 
with a significant barrier separating it from the 
metastable minimum at (j> = 0. The second plot 
displays the profile of the saddle point <f>\,{r) in 
units of kf for the same sequence of scales. For 
k ~ kf the characteristic length scale of the bub- 
ble profile and l/k are comparable. This is an 
indication that we should not proceed to coarse- 
graining scales below kf. 

Our results for the nucleation rate are pre- 
sented in the third plot. The horizontal axis cor- 
responds tok/ sJU k {4>t), i.e. the ratio of the scale 
k to the mass of the field at the absolute mini- 
mum located at (f> t . Typically, when k crosses 
below this mass, the massive fluctuations of the 
field start decoupling. The evolution of the con- 
vex parts of the potential slows down and even- 
tually stops. The dark diamonds give the val- 
ues of the action Sk (free energy rescaled by the 
temperature) of the saddle point at the scale k. 
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We observe a strong k dependence of this quan- 
tity. The stars indicate the values of ln(Ak/kj). 
Again a substantial decrease with decreasing k 
is observed. This is expected, because the fluc- 
tuation determinants in Ak are calculated with 
an effective ultraviolet cutoff of order k. The 
dark squares give our results for — ln(//fci) = 
Sk — \n(Ak/kj). The k dependence of this quan- 
tity almost disappears for kj y/Uj!^) ^ 1- The 
small residual dependence on k can be used in or- 
der to estimate the effect of the next order in the 
expansion around the saddle point. This contri- 
bution is expected to be smaller than \a{Ak/kj). 
It is apparent from fig. 1 that the leading con- 
tribution to the pre-exponential factor increases 
the total nucleation rate. This behaviour, asso- 
ciated with the fluctuations of the field whose 
expectation value serves as the order parameter, 
is observed in multi-field models as well. 

Our results confirm the expectation that the 
nucleation rate should be independent of the scale 
k that we introduced as a calculational tool. They 



also demonstrate that all the configurations of 
the second plot give equivalent descriptions of the 
system, at least for the lower values of k. This 
indicates that the critical bubble should not be 
associated only with the saddle point of the semi- 
classical approximation, whose action is scale de- 
pendent. Only the combination of the saddle 
point and its possible deformations in the ther- 
mal bath has physical meaning. 

For larger values of \k /{~ m fc ) 1 ^ 2 the de- 
pendence of the nucleation rate on k becomes 
more pronounced. We demonstrate this in the 
second row of plots, for which Afc = 0.1 • ko, 
m l = -0-013 • fcg, 7 fco = -1.61 • 1(T 3 kl /2 and 
kf/ko = 0.0421. The dimensionless coupling is 
x k„/(- m l a ) 1/2 = 0-88 (instead of 0.48 for the 
first row). The first-order contribution ln(Afc /kj) 
in the expansion around the saddle point is now 
comparable to the lowest order contribution Sk- 
This indicates that higher orders in the expan- 
sion are important and the series stops converg- 
ing. Therefore, there is a limit to the validity of 
Langer's theory of homogeneous nucleation. 
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3. Validity of the theory of homoge- 
neous nucleation 

An intuitive derivation of the region of appli- 
cability of this theory can be obtained through 
use of approximate analytical expressions for the 
pre-exponential factor. The dark triangles in the 
third plot of each row in fig. 1 display values pre- 
dicted by the approximate expression [4, 7] 

oo -1 1/2 

r 3 [£#(0 6 (r))-E#(O)]dr 

(3.1) 

Good agreement is observed with the numerical 
results denoted by stars, especially in the first 
case in which the expansion around the saddle 
point converges. 

Let us consider a potential given by eq. (2.1), 
with m\ o > 1 . As the bubble- nucleation rate is 
largely independent of k when homogeneous nu- 
cleation theory is applicable, we can obtain an 
estimate of the rate by using k — ko and eq. (3.1). 
Through appropriate rescalings of r and </> we can 
define a dimcnsionless potential U((f>) = (f> 2 /2 — 
3 /3 + with h = 9A feo m| o / 7 | o . For 

ftsil the two minima of the potential have ap- 
proximately equal depth, while for h — > the 
metastable minimum at the origin has negligi- 
ble depth compared to the stable one (spinodal 
point). One can estimate [4] 



R 



^(A k Jk 4 ) 
Sk 



T{h)^. 



(3.2) 



The function T(h) diverges for h — > [4], sig- 
nalling the breakdown of homogeneous nucleation 
theory in the spinodal region. Moreover, this 
breakdown may occur even away from the spin- 
odal region. T(h) is small for h close to 1 [4], but 
a sufficiently large value of \k / m k can lead to 
R <; 1. This case corresponds to weak first-order 
phase transitions, as can be verified by observ- 
ing that the saddle-point action, the location of 
the true vacuum and the difference in free-energy 
density between the minima go to zero in the 
limit mk /\ko —* for fixed h. 

The above conclusions are confirmed by the 
numerical computation of the nucleation rate. 

< are related by 



> and ml 



1 Models with m| 
a field shift [4]. We have employed this shift in order to 
move the metastable minimum to the origin in fig. 1. 



A complete analysis of the region of validity of 
homogeneous nucleation theory for potentials of 
the form of eq. (2.1) was presented in rcf. [4]. 
Also, two-scalar theories with potentials of var- 
ious forms were studied in ref. [5]. The conclu- 
sions can be summarized as follows: Homoge- 
neous nucleation theory relies on an expansion 
around the saddle point that is assumed to dom- 
inate the tunnelling process (the critical bubble). 
Therefore, it is applicable as long as this expan- 
sion converges. Its breakdown is signalled by a 
pre-exponential factor comparable to the leading 
exponential suppression and a strong dependence 
of the predicted nucleation rate on the coarse- 
graining scale. 

4. Radiatively induced first-order 
phase transitions 

These transitions [14] form a class of special in- 
terest, as many cosmological phase transitions 
belong to it. A typical example is presented in 
fig. 2 [5] for a two-scalar theory with an initial 
potential 



9kg 

4 ' 



(4.1) 



with 4>l ko = 1.712 fc , A fco = 0.01 fc and g ka = 
0.2 ko. In the first plot of the first row we depict 
a large part of the evolution of Uk((f>i,<f>2 = 0). 
The initial potential has only one minimum along 
the positive 0i-axis and a maximum at the ori- 
gin. In the sequence of potentials depicted by 
dotted lines we observe the appearance of a new 
minimum at the origin at some point in the evo- 
lution (at kcw/ko ~ e~ 4A ). This minimum is 
generated by the integration of fluctuations of 
the 4>i field, whose mass depends on 4>\ through 
the last term in eq. (4.1) (the Coleman- Weinberg 
mechanism [14]). Immediately below we plot the 
mass term d 2 Uk/ d<p\{<j)i, = 0) of the (f> 2 field 
along the 0i-axis. It turns positive at the origin 
at k = kcw ■ We calculate the nucleation rate us- 
ing the potentials of the last stages of the evolu- 
tion. The solid lines correspond to ki/k = e -4 ' 7 , 
while the line with longest dashes corresponds to 
kf/ko = e- 5 - 2 . 
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Figure 2: A radiatively induced first-order phase transition in a model with initial potential given by eq. (4.1) 
with 4>Q kg = 1.712 ko, Afc = 0.01 fco and gt a = 0.2 fco- The calculation of the bubble-nucleation rate is 
performed between the scales ki = e _4,7 fco and kf = e~ 5 ' 2 ko- All dimensionful quantities are given in units of 
k f . 



In the third plot of the first row we display 
the profile of the saddle point 0&(r) for various 
k between ki and kf. In the middle row we plot 
the quantities 

W lk {r) = ^ (fc = Mr)) ~ |J (<k = 0) 

(4.2) 

for <j>2 = and i = 1,2, i.e. the difference of 
the mass terms of the two fields along the bubble 
profile and the false vacuum. These quantities 
are expressed in units of kf. We observe that 
the mass of the </>2 fluctuations in the interior of 
the critical bubble is much larger than the other 
mass scales of the problem, which are comparable 
to kf. This is a consequence of our choice of 
couplings 3fc /Afc = 20. Such a large ratio is 
necessary for a strong first-order phase transition 
to be radiatively induced. 

The results for the nucleation rate are sum- 
marized in the last plot of the second row. The 
horizontal axis corresponds to the ratio of the 



scale k to the mass of the 4>i field at the absolute 
minimum. The dark diamonds give the negative 
of the action of the saddle point — Su at the scale 
k. The stars indicate the values of the prefactor 
\n(Aik/kj) associated with </>i fluctuations, while 
the triangles denote the prefactor ln^fc of the 
4>2 field. A significant k dependence of all these 
quantities is observed. The dark squares give 
our results for ln(I/kj) = — Sk + \n[AikA2k/kj]. 
This quantity has a strong k dependence, which 
signals the breakdown of the expansion around 
the saddle point. This is expected, as | In A 2 k\ is 
comparable to the saddle-point action Sk, even 
though \n(Aik / k^) remains small. 

The breakdown of the expansion is a generic 
problem of all the radiatively induced first-order 
phase transitions we have studied [5]. It can be 
understood in physical terms by considering the 
quantity W / 2fe(r) defined in eq. (4.2) and depicted 
in the second plot of the second row in fig. 2. Due 
to the large values of the ratio gk /^k = 20, the 
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mass of the fa field is much larger than kf, the 
typical scale of the tunnelling process. This in- 
dicates that fluctuations in the <p2 direction cost 
excessive amounts of free energy and are sup- 
pressed. As these fluctuations are inherent to the 
system, the total nucleation rate should be sup- 
pressed as well. This explains the very negative 
values of In A 2 k ■ This behaviour is not surprising. 
The radiative corrections to the potential and the 
pre-exponential factor have a very similar form 
of fluctuation determinants. When the radiative 
corrections are large enough to modify the ini- 
tial potential and generate a new minimum, the 
pre-exponential factor should be expected to be 
important also. 

We point out that the presence of a region of 
negative values for the function W\k(r), depicted 
in the second plot of the first row in fig. 2, results 
in the existence of a class of <f>\ fluctuations with 
small typical free energy. As a result, deforma- 
tions of the critical bubble in the 4>i direction are 
favourable and their effect tends to enhance the 
total nucleation rate. This explains the positive 
values of \n(Aik/kj). 

Finally, we mention that the prefactor A 2 k 
can be estimated using approximate expressions 
analogous to eq. (3.1) and appropriate rescalings. 
One finds [7] 

(4.3) 

bk f 9k„ 

with R(h) 24 for any h. This expression indi- 
cates that only very small values of Xk /gk may 
lead to a convergent expansion around the sad- 
dle point. However, eq. (4.3) is not valid for Xk 
smaller than gk by more than two orders of mag- 
nitude. We have not found any model for which 
the numerical study of a radiatively induced first- 
order phase transition indicates a convergent ex- 
pansion. 

5. Two dimensions 

A crucial confirmation of the reliability of our ap- 
proach was obtained in ref. [6] , where bubble nu- 
cleation was studied for (2+l)-dimensional theo- 
ries. Our results were compared with those of the 
lattice simulations of ref. [15]. We considered po- 
tentials similar to the one of eq. (2.1), with small 



modifications in order to match the continuum 
theory of the lattice simulations. The technical 
aspects of the application of our formalism to two 
dimensions can be found in ref. [6]. 

In fig. 3 we present a comparison of results 
obtained through our method with lattice results. 
The dimcnsionlcss parameter A is approximately 
equal to 3A/c / (6 m\ ) in terms of the parame- 
ters of the potential of eq. (2.1). The dimension- 
less variable 9 (that plays the role of a rescaled 
temperature) is approximately equal to 7! /mj. . 
The values of 1/9 are given on the horizontal axis 
of the plots in fig. 3. (For the details, see ref. [6].) 
The diamonds denote the saddle-point action Sk ■ 
For every choice of A, 9 we determine Sk at two 
scales: 1.2 kf and 2kf. The light-grey region be- 
tween the corresponding points gives an indica- 
tion of the k dependence Sk- The negative of the 
bubble-nucleation rate — In (//m 3 ) is denoted by 
dark squares. The dark-grey region between the 
values obtained at 1.2 kf and 2k f gives a good 
check of the convergence of the expansion around 
the saddle-point. If this region is thin, the pref- 
actor is in general small and cancels the k de- 
pendence of the action. The dark circles denote 
the results for the nucleation rate from the lattice 
study of ref. [15]. 

For A = 0, 0.1, 0.2 the values of -ln(//m 3 ) 
computed at 1.2 kf and 2 k f are equal to a very 
good approximation. This confirms the conver- 
gence of the expansion around the saddle-point 
and the reliability of the calculation. The k de- 
pendence of the saddle-point action is cancelled 
by the prefactor, so that the total nucleation rate 
is k independent. Moreover, the prefactor is al- 
ways significantly smaller than the saddle-point 
action. The agreement with the lattice predic- 
tions is good. It is clear that the contribution of 
the prefactor is crucial for the correct determina- 
tion of the total bubble-nucleation rate. 

For larger values of A the matching between 
the lattice and the renormalized actions becomes 
imprecise. As a result, we cannot make sure that 
we are studying the same theory as the one simu- 
lated on the lattice. This generates deviations of 
our results from the lattice ones, which start be- 
coming apparent for A = 0.25. However, the in- 
ternal consistency criteria of our method provide 
a test of the convergence of the expansion around 
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Figure 3: Comparison of our method with lattice studies: Diamonds denote the saddle-point action Sk and 
squares the bubble-nucleation rate — In (//m 3 ) for k = 1.2 kf and 2 kf. Dark circles denote the results for the 
nucleation rate from the lattice study of ref. [15]. 



the saddle point in all cases. The consistency of 
our calculation is achieved for 1/9 <^ 0.12 even for 
A = 0.32. However, the breakdown of the expan- 
sion is apparent for 1/9 <, 0.12. The prefactor 
becomes comparable to the saddle-point action 
and the k dependence of the predicted nucleation 
rate is strong. 

6. Applications 

As a first example we consider the bound on the 
Higgs-boson mass from the stability of our vac- 
uum. The top-quark radiative corrections to the 
zero-temperature effective potential of the Higgs 
field may result in the appearance of a new min- 
imum, deeper than the one located at 247 GeV. 
In order for this not to happen, a lower bound 
on the Higgs-boson mass must be imposed. This 
bound can be relaxed if one allows for the pres- 
ence of a new vacuum, but demands that the 
time needed for a transition to it is larger than 
the age of our universe. The largest rates are as- 
sociated with thermal transitions at high temper- 



atures, from a metastable minimum in the sym- 
metric phase of the Standard Model directly to 
the absolute minimum at very large Higgs-field 
expectation values [16]. 

We consider the effective three-dimensional 
description of the top-Higgs system at a scale 
ko ~ T. The top-quark fluctuations are almost 
decoupled because of the absence of a zero Mat- 
subara frequency for fcrmions. This permits a 
very simple determination of the potential Uk (<p), 
which is equal to the sum of the zero-temperature 
potential and the contribution from the top-quark 
thermal fluctuations [7] . It can be written as 

where m 2 = g 2 4 T 2 /4, k = | A4 T, gn is the top- 
quark Yukawa coupling and A4 the average (neg- 
ative) value of the quartic coupling over the re- 
gion of Higgs-field expectation values relevant for 
the critical-bubble profile. The term ~ (j) 2 is the 
leading temperature correction in the expansion 
of the one-loop contribution from the top quark. 
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The action of the saddle point is [7] 

S k0 ~ 18.9-= 9.45 i" (6.2) 
k \M\ 

The prefactor can be estimated through eq. (3.2) 

ln(^ o /fc 4 )«2.00^ = ^. (6.3) 

For the experimental value of the top-quark mass 
{gu ~ 1) we nn d | -X4 1 ^ 0.05 for all Higgs-field 
expectation values, in agreement with ref. [16]. 
As a result, the prefactor is expected to give only 
a small correction to the bubble-nucleation rate. 

The electroweak phase transition is radia- 
tivcly induced by gauge field fluctuations. Even 
though the generalization of our formalism to 
gauged systems requires further work, an esti- 
mate of the validity of the expansion around the 
saddle point can be obtained by exploiting the 
resemblance between the gauge-field fluctuation 
determinants and the ones in two-scalar models. 
The estimate for the prefactor \\nA 2 k f \, associ- 
ated with the gauge fields, is given by eq. (4.3), 
with X/g taking the "effective" value [7] 

/A\ ^ (4m w + 2m z ) m 2 H = / m H \ 2 

Wsm~ 4(4m^ + 2m|) ' VlOOGeW 

(6.4) 

For the prefactor to be smaller than 1/2 of the 
action, m# Ss 25 GeV is required. 

Higgs-boson masses below the experimental 
lower limit of 90 GeV are of academic interest. 
Moreover, for tuh <5 80 GeV, there is no phase 
transition in the high-temperature Standard Mo- 
del. There is the possibility, however, for a suffi- 
ciently strong first-order phase transition for ba- 
ryogenesis within the Minimal Supcrsymmctric 
Standard Model with a very light stop [17]. In 
this model, the "effective" value of g/X obeys [7] 

/A\ (Am w + 2m z + 6m t ) m 2 H 

V 9 /mssm 4(4to=V + 2m% + 6m?) 

= »■" (To^v) 2 - <«> 

For m t = 175 GeV, in order to have a prefactor 
smaller than 1/2 of the action, mjj ^ 45 GeV is 
required. 
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